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Imperfect photometric calibration of galaxy surveys due to either astrophysical or instrumen- 
tal effects leads to biases in measuring galaxy clustering and the resulting cosmological parameter 
measurements. More interestingly (and disturbingly), the spatially varying calibration also gener- 
ically leads to violations of statistical isotropy of the galaxy clustering signal. Here we develop, 
for the first time, a formalism to propagate the effects of photometric calibration variations with 
arbitrary spatial dependence across the sky to the observed power spectra and to the cosmological 
parameter constraints. We develop an end-to-end pipeline to study the effects of calibration, and il- 
lustrate our results using specific examples including Galactic dust extinction and survey-dependent 
magnitude limits as a function of zenith angle of the telescope. We establish requirements on the 
control of calibration so that it doesn't significantly bias constraints on dark energy and primordial 
non-Gaussianity. Two principal findings are: 1) largest-angle photometric calibration variations 
(dipole, quadrupole and a few more modes, though not the monopole) are the most damaging, and 
2) calibration will need to be understood at the ~ 0.1%-1% level (i.e. rms variations mapped out 
to accuracy between 0.001 and 0.01 mag), though the precise requirement strongly depends on the 
faint-end slope of the luminosity function and the redshift distribution of galaxies in the survey. 



I. INTRODUCTION 

Large-scale structure (LSS) measurements have be- 
come an extremely powerful probe of cosmology over the 
past 30 years. Starting with the pioneering Harvard-CfA 
survey [9], all the way to the Sloan Digital Sky Survey 
pH] and its extension Baryon Oscillation Sky Survey [5], 
Two-degree Field survey [5], and WiggleZ [IT] , the LSS 
surveys have revolutionized our understanding of the dis- 
tribution of matter and energy in the cosmos, and helped 
impose percent-level constraints on the cosmological pa- 
rameters (e.g. [3]). 

A major challenge in current and future LSS surveys 
is understanding of the sample selection. We define cal- 
ibration to be the measure of our understanding of the 
selection of our sample of galaxies, and calibration errors 
to be any unaccounted-for angular and redshift variations 
in the selection. The purpose of this paper is to deter- 
mine how well calibration errors need to be controlled in 
order to avoid substantial degradation of the information 
we can extract from the LSS. 

A particular source of uncertainty is known as photo- 
metric calibration. The term refers to the adjustments 
required to establish a consistent spatial and temporal 
measurement of flux of the target objects in the different 
bands of observation throughout the entire photometric 
survey. This is an enormous problem that all existing and 
upcoming wide area surveys face. The difficulty comes 
from the variability of various building blocks of the ob- 
servational pipeline, which makes it difficult to establish 
a consistent flux baseline at each band (i.e. the flux ze- 
ropoints). In other words, because the instrument sensi- 
tivity is constantly changing, and so are the sources and 
intensity of noise, it is difficult to consistently compare 



the fluxes for objects at different parts of the sky imaged 
at different times. Some examples of the manifestations 
of the photometric calibration errors in surveys are: 

• Detector sensitivity: At any given time, the sensi- 
tivity of the pixels on the camera vary along the 
focal plane. In addition, the sensitivity of a given 
pixel can change with time. 

• Observing conditions: Wide area surveys are car- 
ried out over several years and conditions are con- 
stantly changing. Observing conditions suffer from 
constant spatial and temporal variations. 

• Bright objects: The light from foreground bright 
stars and galaxies affects the sky subtraction proce- 
dure, which impairs the surveys' completeness near 
bright objects [U [28], and distorts the measured 
shapes of these faint galaxies [55] . 

• Dust extinction: Dust in the Milky Way absorbs 
light from the distant galaxies. As we show later, 
imperfect extinction correction can have serious 
consequences to cosmological clustering analyses. 

• Star-galaxy separation: In photometric surveys, 
faint stars can be erroneously included in the galaxy 
sample. Conversely, galaxies are sometimes mis- 
classified as stars and culled from the sample. 
These effects are important because stars are not 
randomly distributed across the sky. 

• Deblending: Galaxy images can overlap, and it can 
be difficult to cleanly separate photometric and 
spectroscopic measurements for the blended ob- 
jects. 
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Variabilities in the instrument sensitivity and observing 
conditions cause an angular variability in the depth of 
observations that the survey can achieve through each 
filter. Variations in the depth result in angular varia- 
tions in the number density and redshift distribution of 
objects. In addition, because galaxy spectra are not flat, 
and because the sample selection involves more than one 
filter, depth variations cause variations in the angular 
and redshift distribution of galaxy types. 

This variability in the sample selection can, in princi- 
ple, be accounted for. This is not always done however, 
and it is common, for example, for correlation analyses 
of current data to assume a constant depth for the entire 
survey. Indeed, several sources of variability have been 
accounted for in the analysis of existing data - see in par- 
ticular the pioneering work on the subject in the modern 
era of LSS surveys by Scranton et al. [31] ( see also Voge- 
ley [33]), and the more recent efforts by Ho et al. [17] and 
Ross et al. [28] . These authors modeled a wide variety of 
systematic errors, some of which qualify as the calibra- 
tion errors (e.g. seeing, airmass, calibration offsets). In 
particular, the latter two papers identified bright stars as 
the major contaminant which adds significant power to 
the intrinsic clustering signal at large scales, and they ap- 
plied two separate successful techniques to subtract this 
systematic contamination. 

For the upcoming surveys an even more detailed anal- 
ysis will be needed, ideally utilizing a formalism that 
is suited to a wide variety of photometric calibration 
systematics mentioned above and captures any kind of 
calibration-related systematic. One would also like to 
provide a guidance on how much calibration error, as a 
function of scale, can be tolerated in order not to degrade 
the cosmological parameter inferences. Here we aim to 
address both of these desiderata. 



II. GOALS OF THIS WORK 

In this paper we set out to study calibration errors in 
the most general way possible. Our goal is to build an 
end-to-end pipeline into which we can feed in calibra- 
tion errors (or uncertainties) due to an arbitrary cause, 
and from which we obtain biases in cosmological param- 
eters inferred from measurements of galaxy clustering in 
some LSS survey. We then turn the problem around, and 
estimate how well the calibration errors need to be con- 
trolled in order not to appreciably bias the cosmological 
parameter estimates. 

To keep the scope of this paper reasonable, we only 
consider measurements of the galaxy two-point correla- 
tion function (i.e. the power spectrum), and leave other 
observable quantities - higher-order correlation functions 
of galaxies, for example - for future work. We also do 
not consider the effect of the photometric redshift errors 
which, while very important, are not expected to change 
our results in a major way, so we leave the photo-zs for 
a future analysis. 



The paper is organized as follows, in Sec. 



Ill we de- 
under- 



scribe our formalism of modeling both the true, 
lying galaxy density field and the systematic errors de- 
scribing variations in the photometric calibration. In 



Sec. IV we present the formalism to derive cosmological 
constraints and biases on cosmological parameters. In 
Sec. [V] we propagate the effects of the systematic errors 
to calculate the biases in the cosmological parameters. 
We conclude in Sec. |Vl| Important but technical details 



regarding various aspects of the computation of the ef- 
fects of the photometric variation systematics on the ob- 
servable quantities are relegated to the four Appendices. 



III. FORMALISM: DESCRIBING SPATIALLY 
VARYING CALIBRATION 

In this section we start by defining calibration errors 
and their field c(n), and proceed to derive the biased 
galaxy fluctuations in terms of this field in multipolc 
space. 



A. Calibration errors: definition and basics 

Let true galaxy counts on the sky be denoted by N(ii), 
where n is an arbitrary spatial direction. The survey 
mean is given by N = (iV(n)) s k y , where the average here 
is taken over the observed sky. These true fluctuations 
in the galaxy counts can be expanded into harmonic co- 
efficients ag m as 



N(h) - N 
N 



EE 

i=0 m =-l 



(1) 



Consider a survey where a deterministic calibration er- 
ror c(n) biases galaxy counts. In other words, given the 
true galaxy number counts in some direction N(n), the 
observed number is 



JV obB (n) = [l + c(n)]iV(n), 



(2) 



which implicitly defines the calibration field c(n). We can 
expand the calibration field relative to its fiducial value 
of zero (corresponding to no error) 



c(n) 



£ 1= m x = -k 



Qimi ^imi 



(3) 



where hereafter we assume that the calibration error 
dominates on large scales, and persists only out to some 
maximum multipole ^ ca iib,max, corresponding to angular 
scale of 7r/4aiib,max radians. 

The statistical properties of the two galaxy number- 
density field, and the calibration-error field are, respec- 
tively 



[Cir, 



0: 



{at m a\, m t) — S mm '5wCg (4) 
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Throughout the paper, angular brackets (•) indicate en- 
semble averages, that is, averages over different realiza- 
tions of the Universe. To reiterate, N(h) is the Gaussian 
random, isotropic field as predicted by inflation, while 
c(n) is a deterministic function given by calibration er- 
rors in the survey at hand. 

In the remainder of this paper, we use the following 
definition: Calibration variations (or errors) are depar- 
tures o/c(n), or its harmonic coefficients ci m , from zero. 
Our goal is to estimate how accurately those variations 
have to be known in order not to bias the cosmological 
parameter estimates. 

Notice that we do not lose any generality by assuming 
that the c-ficld is fixed, rather than stochastic like the 
true galaxy density field. In this paper, we are effectively 
asking how much does this fixed systematic error bias 
the usual cosmological constraints. We are, of course, 
free to iterate over a number of specific incarnations of 
this 'fixed' error. A more specific example would be to 
ask how much does the Galactic dust pattern - its di- 
rection and amplitude fixed for the moment - bias some 
cosmological inference if unaccounted for perfectly, and 
then to repeat the analysis for a number of dust pattern 
realizations, or even different dust models. 



B. Galaxy clustering and calibration errors: 
general case 

In this Section we present the main results regarding 
the effect of the calibration errors on the observed clus- 
tering of galaxies and leave the details of the derivation 
to Appendix [A] 

The observed (i.e. biased) mean and density fluctua- 
tion are given by (see Eq. (|A6|): 



where 



e= ^(c(n)7V(n)) sky . 



(7) 



The observed galaxy overdensity field can also be ex- 
panded in terms of the harmonic basis 



,5 obs (fi) = t(fi) = ^W^ m (n) 



(8) 



tin. 



The harmonic coefficients of S ohs are calculated to be (see 
again Appendix fAl 



tin 



1 



f 



at,. 



%m ! m c <imi B < 2 m 2 



(9) 



where we dropped the I — 0, to — term (see Eq. (AfO) 
for the full expression) since we can not use the galaxy- 
fluctuation monopole to determine the cosmological pa- 
rameters, The coupling matrix R can be written in terms 
of Wigner 3j symbols 



= (-iy 



(2*i + l)(2*a + l)(2* + l) 



47T 



h 4 A ( l\ I 




mi nri2 —to 



(10) 



6° hs (h) = — [«S(n)(l + c (n)) + c(n)- 
N ohs = N(l + e) 



The ensemble average of the multipole moments be- 
(6) comes, after some algebra 



(1 



5w Ci + 



isotropic 



breaks statistical isotropy 



(11) 



where we defined 

U i3 



Qimi R 



(12) 



which is a function that depends on the Wigner 3j sym- 
bols as well as the calibration- field coefficients ct m . 



Equation (11) is the key result in this paper. As the 



label in the equation shows, the observed galaxy density 
field t(h) exhibits broken statistical isotropy. In partic- 
ular, the variance of t is not rotationally invariant any 
more (i.e. it depends on to), and covariance between the 
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different £ modes is not zero any more. We can now uti- 
fize this formula and consider the isotropically measured 
power (i.e. assuming £ = i! and averaging over to = m!) 
and estimate how accurately any given systematic, de- 
scribed by the full set of q,„, needs to be understood 
in order not to degrade the accuracy in measuring the 
cosmological parameters including non-Gaussianity 1 . 



C. Galaxy clustering and calibration errors: 
isotropic power case 

We usually - essentially always, in fact! - assume that 
the field is isotropic, and then we use the data to calculate 
the correlation function, power spectrum, etc. Let us 
see how the assumed-isotropic angular power spectrum 



is biased in terms of an arbitrary contamination field. 



Setting £ = £' and m = m' in Eq. (11), we get 



(|W| 2 > — — ^ 2 (Smm'See'Ce + 2 (l/" m ) C e + 

+ Ei^r^+w 2 ) as) 



To assume statistical isotropy, we not only set £ = £' and 
m = m' but further average over the 21+1 values of m 
for a fixed I. Then we obtain the prediction for the an- 
gular power spectrum that one would measure assuming 
statistical isotropy even when the systematics break it: 



Ti 



2£+l 



coo 



(1 + ^) 2 



471\ 



i •• \T7 e '~ 

= — (. \Z—i£ 2 m2 I msm 



Cf„ 



2£+l 



(14) 



where \U\ 2 = UU* , and where the sum over £2 goes in 
principle over all multipoles (though only those from the 
range [•£-4alib,max^ + 4alib,mad are nonzero), while to 2 
goes from —£2 to £2- Note that the term linear in U seen 



Tt 



in Eq. ( |13| dramatically simplified in the expression for 
(Eq. 



14 ) after we used the summation relation 



£ 

m=-t 



(-l) m ( £l 1 M = (-l)'V2£+T^ l0 . (15) 
v ' V mi to —to I 1 v ' 



For a pure monopole calibration error (i.e. a pure coo 
term) , one can verify that the effects of the e term and the 
cqo term in Eq. ( 14 ) exactly cancel and is unchanged. 



This makes intuitive sense, as a shift in the monopole 
changes the mean counts on the sky but does not affect 
the density fluctuations. 

One can intuitively understand terms in Eq. ( 14 ) : 



The (1 + e)~ 2 prefactor accounts for the change 
in the mean number of galaxies observed on the 
sky; see Eq. (A7). As already mentioned, this term 



effectively ensures that the effect of the monopole 
(cqo, or constant change in calibration across the 
sky) is precisely canceled out. 

• The terms containing U^ m introduce coupling be- 
tween the multipoles. In particular, if there are 
calibration error multipoles out to some multipolc 



^caiib,max, then there will be contributions in the 
range [i-4alib,max, ^+4alib,max] at a given I, a fact 
familiar from spin coupling in quantum mechanics. 

• Presence of the term |Q m | 2 essentially means that 
the power of the calibration field is added to the 
intrinsic galaxy power spectrum. In other words, 
even if the distribution of galaxies on the sky were 
perfectly uniform so that 5 = N/N —1 = over 
some area, the calibration field will induce power 
so that S ohs ^ 0. 



From the structure of Eq. ( 14 ) , it is clear that calculat 



ing and storing the coefficients U is challenging. Naively, 
the problem requires calculation of the order of 10 18 co- 
efficients. Appendix [D] describes our approach of lim- 
iting the number of evaluations of (£, t\, £2, to, to 1? 777.2) 
and tabulating the coefficients U so that the number of 
operations is only of order 10 8 (for £ max = 1000 binned 
in ~ 30 multipole bins and considering the calibration 
variations out to 4aiib,max = 20), and is thus feasible. 

We plot the biased power spectra Tg further below in 
the next section. 

Finally, it is worth writing down the observed angular 
cross-correlation power spectrum between fluctuations 
t)^ and t*^ in two different tomographic redshift bins i 



and j; it follows straightforwardly from Eq. (14 1 that 
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(ij) 



(i) 



2€4 



(16) 



where c^, and (U^ m )^ are all evaluated in the 

redshift bin z (and same for j), and where are the 

cross-correlation power spectra. While physical sources 
of calibration error are typically local and thus redshift - 
independent, in Sec. [V] we demonstrate that converting 
from the magnitude error to the calibration field c(n) = 
(SN/N)(n) depends on the faint-end slope of the lumi- 
nosity function, which typically is redshift-dependent, 
hence making the harmonic coefficients of c(n) also z- 
dependent and thereby potentially introducing couplings 
between the different redshift bins. 



D. Additive and multiplicative systematics 

Before we obtain numerical results on how some real- 
istic calibration errors affect the observed power spectra, 
it pays to consider qualitatively how the power spectrum 
of galaxies is affected. 

It is often useful to divide the effect of systematic errors 
into additive (those whose field is added to the true field 
observed on the sky), and multiplicative (those whose 
field multiplies the true field); see e.g. [13] and [15] |2"2"] 
for the CMB and weak lensing systematic examples, re- 
spectively, where this nomenclature has been previously 
employed. These terms refer to the systematic error that 
either adds to the true galaxy fluctuations, or else mul- 
tiplies it and modulates the true signal; see Eq. ([6| for 
the real-space and Eq. ([9| for the harmonic-space pic- 
ture. For example, the angular parentheses on the right- 
hand side of Eq. ^ contain three terms; the term ag m 
corresponds to the true density field, C£ m represents the 
additive effect of the systematic error, while the term con- 
taining ci m ai m term together with the geometric factor 
R and the appropriate sum, represents the multiplicative 
effect of the systematics 2 . 

Additive and multiplicative errors in the counts trans- 
late into additive and multiplicative contributions from 
the calibration field to the observed galaxy power spec- 
trum, see Eq. ( 11 ): the additive error is the term ct m c* tm 



while the multiplicative error are all terms involving the 
coupling matrix U . The two kinds of errors produce qual- 



itatively different effects: additive error at some multi- 
pole £\ only affects power at that multipole, while the 
multiplicative error affects power at a range of multi- 
poles; in particular, true power at an arbitrary multipole 
£ would leak to all multipoles in the range [£ — £%, £ + £%]. 

The additive terms dominate the error budget on the 
largest scales, but are subdominant at smaller scales [33] , 
This can be understood qualitatively as follows: modulo 
geometric coupling terms, both additive and multiplica- 
tive terms are proportional to the square of coefficients 
C£ m , but the multiplicative terms are further multiplied 
by the fiducial angular power Ci. Given that C( « 1 
at all £ and any redshift 3 , the multiplicative terms are 
suppressed relative to the additive terms. At higher mul- 
tipoles, on the other hand, there are more ways in which 
power from other scales can leak into that I so that the 
sums associated with multiplicative terms make them the 
dominant systematic contribution. 

While it is often assumed for simplicity that the sys- 
tematic errors, calibration or other, are represented by 
purely additive errors, we just demonstrated that both 
kinds of errors are important. In fact, in the plausible 
scenario where largest-scale information in the survey is 
ignored to avoid the systematic contamination, the mul- 
tiplicative errors dominate. In what follows we use the 
full expressions containing both additive and multiplica- 
tive terms. 



IV. POWER SPECTRA AND THEIR BIASES 

In this section we propagate the effect of calibration 
errors to estimate biases in the cosmological parameters 
describing dark energy and primordial non-Gaussianity. 



A. Fiducial cosmological model 

We consider a set of cosmological parameters with the 
following fiducial values: matter density relative to criti- 
cal f2jvf = 0.25, equation of state parameter today wq = 



2 Technically speaking, the systematic effects are all multiplicative 
in Eq. (JsJ because of the 1/(1 + e) prefactor; however given that 
this prefactor is typically close to unity, the second and third 
term in parentheses of this equation act approximately in addi- 
tive and multiplicative sense. 



3 One exception - C{ of order unity or larger — is realized in the 
scenario of galaxies at a very low redshift, z < 0.05. However 
such a sample would probably not be useful for cosmology given 
the significant galaxy peculiar velocities and the necessarily small 
volume probed. Moreover, the calibration errors would have less 
impact to begin with, since they would be affecting a very large 
intrinsic clustering signal. 
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— 1, its variation with scale factor w a = 0, spectral index 
n = 0.96, and amplitude of the matter power spectrum 
In A where A = 2.3 x 10~ 9 (corresponding to a s = 0.80) 
defined at scale scale k = 0.002Mpc~ 1 . Note that we hold 
fixed the Hubble constant h — ffo/(100km/s/Mpc) (or 
equivalently, physical matter density Qm^ 2 ), the physi- 
cal baryon density Qgh 2 , and we assumed a flat universe. 
On the other hand, we do not assume any other prior in- 
formation, such as the CMB information from WMAP 
and Planck. In practice, this prior information would 
largely serve to fix the aforementioned three parameters. 
Note that this rather restricted set of assumptions about 
the set of cosmological parameters and external infor- 
mation about them is sufficient for our analysis: we are 
primarily concerned about the effect of the calibration 
systematics on the measured power spectra, and on the 
biases in the dark energy and non-Gaussianity parame- 
ters. Given that the systematics strongly depend on the 
properties of the galaxy sample and the survey (as we 
discuss further below), it is not necessary to model the 
up-to-date knowledge about the cosmological parameters 
in great detail. 

We assume a survey covering 5,000 square degrees (so 
/sky — 0.12) with information out to z max = 1, cor- 
responding roughly to the Dark Energy Survey (DES). 
We assume that the number density of the galaxies is 
n(z) oc z 2 exp(— z/zq) with zo = 0.3, and that photo- 
metric redshifts enable splitting the sample in five tomo- 
graphic bins centered at z = 0.1, 0.3, 0.5, 0.7 and 0.9; 
see the left panel of Fig. [T] The fiducial statistical con- 
straints on the dark energy parameters are o~(wo) = 0.06 
and <r(w a ) = 0.24. 

Instead of the power spectrum in wavenumber P(k), we 
consider measurements of the angular power spectrum of 
galaxy fluctuations. In the Limber approximation, which 
is valid on intermediate to small angular scales, the an- 
gular power is given as (e.g. [2"Tl [55] ) 



P. 



(W) — /„(») *0') 



\ a lm a tm I 



(17) 



2tt 



r(z)H(z)b 2 (k, z)A 2 {k, z)W i (z)W j (z)dz 



where i and j are referring to one of the five redshift 
bins 4 . The full, beyond-Limber expression, as well as the 



4 While we use the all of the cross-correlations, i > j, we notice 
that the cosmological constraints - and requirements on the sys- 
tematics - remain unchanged if we only use the angular power 
spectrum auto-correlations with i = j. This is not too surprising, 
as auto-correlations C) are guaranteed to be zero in the Limber 
approximation and in the absence of photometric redshift errors 
(photo-zs generically lead to overlap of redshift bins and thus bin- 
to- bin correlations). At large scales, where we do not employ the 
Limber approximation where it is not accurate (especially for 
nonzero /nl which introduces very large-scale correlations), we 
explicitly verify that cross-powers do not appreciably alter our 
requirements on the control of the calibration errors. 



definition of the W(z) terms, are given in Appendix |b| 
The power spectrum A 2 (fc, z) = k 3 P(k)/ (2ir 2 ) is calcu- 
lated using the transfer function output by CAMB, and 
its nonlinearitics are modeled with the Smith et al. [32 
formulae that were based on a halo model and fit to sim- 
ulations. Appendix [B] has all of the details of how we 
calculate the power spectrum. 

We consider information from multipoles 1 < £ < 1000, 
corresponding to spatial scales from about 10 arcmin 
to 180 degrees. To obtain accurate constraints on the 
parameter /nl which come from large angular scales, 
we use every individual multipole between £ = 1 and 
^caiib.max = 20; beyond this we use ten more widely sep- 
arated bins with A£ ~ 100. Therefore, we use a total 
of 30 bins in £; for the low-^ ones we do not assume 
the Limber approximation and use Eqs. (B4|, while for 



the higher multipoles we use the Limber approximation, 
Eq. (fl7|) above. 



Finally, we also allow for the presence of primordial 
non-Gaussianity. We adopt the widely studied 'local' 
model of non-Gaussianity 



$(x) = <D G (x) + / NL ($ 2 ; -(<i> 2 ,)), 



(18) 



(where $ is the primordial Newtonian gravitational po- 
tential, <5>g is its Gaussian component, and /nl is a di- 
mensionless parameter), the bias becomes scale depen- 
dent, with a new term that goes as k~ 2 [7] 



b(k) = b a + h L (bo - l)S c 



3ft M Hq 



a g(a)T(k)c 2 k 2 ' 



(19) 



where b is the usual Gaussian bias (on large scales, where 
it is constant), 5 C « 1.686 is the collapse threshold, a is 
the scale factor, £Im is the matter density relative to 
critical, Hq is the Hubble constant, k is the wavenumber, 
T(k) is the transfer function, g(a) is the growth suppres- 
sion factor, and c is the speed of light 5 . We assume the 
fiducial model with the Gaussian bias bo = 2 and zero 
non-Gaussianity, /nl = 0. 

The full set of cosmological parameters that we use is 
therefore 



Pa € {Cl M ,Wo,w a ,n,A,f N i,}. 



(20) 



Given a systematic error in the power-spectrum, 8Ci 1 
the resulting bias in the cosmological parameters, 8p a 
can be estimated using the Fisher matrix formalism as 
follows: 



5p a = 



j 



F, 



ah 



dd a) Cov" 



dC, 



(/?) 



d Pj 



(21) 



This formula has higher-order corrections calibrated to N-body 
simulations and also derivable from theory, but these are typi- 
cally small and not crucial for the present analysis so we ignore 
them here. 
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FIG. 1: Left panel: density distribution of galaxies assumed in this paper, and boundaries of the five redshift bins. We ignore 
information at z > 1, thus roughly modeling the difficulties with establishing accurate photometric redshifts at that range (for 
the DES). Right panel: Angular power spectra for five redshift bins [Cross-correlations between the bins, while used in the 
analysis, are very small and not important nor shown in the figure.]. For the first and fifth bin we show, at the low multipole 
end, the full expression that we use at large scales (see Eq. (B4l) and, at the high-multipole end, the linear power spectrum 
for reference. For the first and fifth redshift bin we also show the cosmic variance errors plus shot noise. 



where a and /? stand for all pairs of bin indices 
with i < j (since C^' = Cf^) and the observed power 
spectrum is equal to the raw power plus shot noise 



(ij) 



P, 



(U) 



Si 



Nf 



(22) 



where A^ sr is the number of galaxies per steradian in the 
tomographic bin i. Moreover, SCi is calculated below as 
the difference between the observed, biased power spec- 
trum Tg and the fiducial CV Assuming observations in 
the linear (and therefore Gaussian) regime only, the co- 
variance matrix is the standard expression that follows 
using Wick's theorem: 



Cov 



r>W) r<( fe 



[21 + 1) / sky A£ 



(23) 



Using these equations, one can easily calculate the bi- 
ases in the cosmological parameters for any given bias in 
the observed angular power spectra. 



B. Biases in the observed power spectra 

We would like to fairly compare the biases as a function 
of £\, so we choose to adopt coefficients describing the 
uncertainties in the calibration field ci imi that lead to 
a fixed variance in the calibration pattern on the sky 
c(n) (and corresponding to, as we will shortly see, fixed 
variance in the angular variations of the magnitude limits 



of the survey). Thus, we have 
Var(c(n)) = 



2ll + 1 ^sys _ |ciimi| 



4tt 



4tt 



(no sum) (24) 



since Cf^ s = \cg imi | 2 / + 1) is the angular power spec- 
trum of the systematics (and really just the sum of their 
coefficients squared) and where the reader is reminded 
that, in this particular calculation we are "turning on" 
one mi) pair at a time. 

To consider the calibration variation in a single multi- 
pole (£i, mi) that leads to a fixed variance in c(n), we 
therefore make the choice 



Rc.Im 



\/47rVar(c(n)) (m = 0) 
V27rVar(c(n)) (m ^ 0) 



(25) 



where in the m ^ case both real and imaginary part of 
the cg imi have the given value. 

Figure [2] shows the difference between the observed 
isotropic part of the power spectrum (see Eq. (14)) 
and the fiducial Ci, divided by the statistical error (cos- 
mic variance plus shot noise) for our assumed DES-type 
survey with / s k y — 1/8. We assume a constant calibra- 
tion error with rms 6 of 0.01 or 0.001 per each multipole 
li separately. While a fixed t\ of the systematic errors 
affects all multipoles £, it affects £ — £i the most, so in 
this graph we only plot the effect on the observed power 



For the purposes of making FigureJSJwe put all of the photometric 
calibration variation in a single value of mi = 0. 
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spectrum at the same multipole at which the systematic 
errors occurs. In other words, Figure [2] shows the max- 
imally affected multipole £, for a fixed calibration varia- 
tion error. The error in the measured tomographic an- 
gular power spectra decreases with multipole £, which 
can be understood easily as follows: to a good approx- 
imation, the additional terms in the observed isotropic 



power spectrum, Eq. ( 14 ) , are dominated by the additive 



term \ci m \^ / (21+ 1). We find this is the case even for the 
lowest-redshift tomographic bin where Ct is the highest 
and thus helps boost the multiplicative terms that con- 
tain the U coefficients. Moreover, generally we find that 
e ~ 0. Therefore, 



Qv, 



Tt — Ct 
<r(C t ) 



(2£ + l) 



{21 + l)/ sky 



Ct 



1 



(26) 



(ii) 

where Ct = C\ is referring to the power spectrum in 
some redshift bin i, and same for and Q m . For a fixed 
variance in the photometric variations, cg m is indepen- 
dent of £, so that this expression goes as 1/ \\j2i + lCg] , 
decreasing with I very approximately as £ -1 ' 3 , at least 
out to £ = 10 plotted in this Figure. 

Figure [2] further shows that higher redshift bins are 
affected more than the lower redshifts. This is easy to 
understand: at higher redshift, the cosmological signal 
is smaller, given that it averages over more large-scale 
structure along the line of sight, and therefore it is more 
susceptible to the (redshift-independent) calibration bias. 

The one important thing to take away from Fig. [2] 
is therefore that the calibration error is expected to 
be most damaging to the deepest surveys (or highest- 
redshift slices of a survey) . And it is precisely those high- 
est redshifts that are most valuable in their information 
about dark energy and primordial non-Gaussianity. 



V. SENSITIVITY TO CALIBRATION ERRORS 

Let us now consider a few specific examples. First, 
we will study the sensitivities to an arbitrary system- 
atic bias in calibration at each multipole separately, i.e. 
one (£i, mi) pair at a time in the Q lOTl . Then we study 
two concrete examples of effects that cause biases in cal- 
ibration: corrections to the dust extinction maps, and 
variable survey depth. 

We first illustrate the effect that a constant-magnitude 
calibration error with some spatial scale has on the biases 
in cosmological parameters. We then give two concrete 
examples of how the magnitude-limit variations can oc- 
cur, and estimate their effect on the cosmological param- 
eters using the formalism that we developed. 
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FIG. 2: Difference between the the observed isotropic part of 
the power spectrum Tt (see Eq. ( |14| ) and the fiducial Ci, di- 
vided by the statistical error (cosmic variance plus shot noise) 
for our assumed DES-type survey with / s k y — 1/8. We as- 
sume a constant rms photometric variation of 0.01 (dashed 
curves) or 0.001 (solid curves) per each multipole I. Note 
that higher redshift bins are affected more than the lower 
redshifts. The fall-off with I can be understood analytically; 
see text for details. 



A. From calibration to galaxy counts 

Consider the observed angular density of galaxies in 

some direction in the sky N(ri) = J n(z, h)dz, where 

n(z, n) is the galaxy density in that direction and at red- 
shift z. 

Calibration errors correspond to variations in the mag- 
nitude limit of the survey (5m max (n). The observed den- 
sity of galaxies changes since the galaxy density is a 
strong function of the survey depth (i.e. the magnitude 
limit). Suppressing the direction label n, we have 



S[log 10 N(z, > m)] 



dlog 10 N{z, > m) 
dm 

s(z)6m niax 



(27) 



where m max is the maximal apparent magnitude observed 
in that direction in some waveband. It follows that the 
systematic bias in the observed fluctuations is 



6N \ 



ln(10)s(z) 5m„ 



(28) 



Often we have information about the selective ex- 
tinction Eb-v\ t ne relation to magnitude extinction is 
5m = 8 A = S[REb~v] where R is the ratio of total to se- 
lective extinction and A is the alternative notation some- 
times used for extinction. Assuming that R is known 
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■ magnitude rms error = 0.01 

(dashed) 
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w=const 




0.001 

(solid) 



4 6 8 10 12 
Multipole of systematics 1 : 



FIG. 3: Bias divided by marginalized statistical error in the cosmological parameters for the fixed magnitude root-mean-squared 
variation of 0.001 (solid curves) or 0.01 (dashed curves), as a function of multipole at which the systematics is introduced. To 
convert the magnitude variation to the SN/N error, we used Eq. ( |28[ ) and the best-fit faint-end slope of the luminosity function, 
s(z), estimated from simulations of [23]; see text for details. The error bars show dependence on which of the m-values contains 
the calibration error; this dependence is small at the largest scales where the calibration error clearly has the largest effect. As 
discussed in the text, the monopole l\ = has no effect on the biases by definition. The line labeled FoM 1 ^ 2 refers to the square 
root of the DETF figure of merit, and serves to gauge any additional dependence brought forth by the temporal variation in 
the equation of state w a . The dashed horizontal line denotes a fixed bias/error ratio of 0.3, which is approximately the upper 
limit of how much effect a systematic error should have on the cosmological parameters without seriously affecting the overall 
constraints in a survey. 



perfectly 7 , 5m ~ RS(Eb-v), an( l thus (restoring the di- 
rection n explicitly) 



^ J (n) = c(n) = ln(10)s(z)5m max (n) (29) 

/ sys 



]n(10)s(z)R6(E B _ v )(n).(30) 



While s(z) is galaxy-population dependent, we can 
still estimate (8N/N) sys to be very roughly of order 
5(Eb-v)(p), given that s(z) is of order 0.1-1 while R 
takes values between about 1 and 5 depending on the 
band; see e.g. the appendix of Schlegel et al. [3U] and 
important updates given in the Table 6 of Schlany and 
Finkbeiner 1251. 



7 If il is not perfectly known, as is often the case, then the mag- 
nitude variation is equal to the variation in the product between 
R and Eg_y. 



B. Calibration Bias per Multipole 

Let us consider biases in our six cosmological param- 
eters as a function of bias in a single multipole C£ imi . 
Following the prescription in the previous section, we as- 
sume that the variance of the calibration field is fixed and 
constant separately at each multipole l\\ see Eq. (l25|). 



We now propagate the calibration variation in a given 
(£i,mi), separately for 1 < l\ < 20 and —l\ <mi < l\, 
and for magnitude given in the above equations, to the 
observed angular power spectra via Eq. (14). 



Figure [3] shows the bias divided by the statistical error 
in cosmological parameters for the fixed magnitude rms 
variation per multipole of (Sm^^tly = 0-01 or 0.001. 
We use Eq. ( 29 1 to translate this to the calibration-field 



variation, and then a modified version of Eq. (25) to cal- 
culate the harmonic coefficients C£ imi that enter the cal- 
culation: 

V47rVar(c(n))/(24 + 1) (m = 0) 

^Re,Im I (31) 

' \/27rVar(c(n))/(24 + 1) (to ^ 0) 
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Maximum multipole of systematics 1, 

r J 1 , max 



FIG. 4: Same as Fig. [3] except now the fixed magnitude error of 0.001 is shared equally among all multipoles in the range 
1 < £i < ^i,max. Under this more realistic assumption than the fixed-multipole of error case shown in Fig. [3] the biases appear 
larger because the contributions from the largest scales dominate the budget at a fixed £i, ma x- 



where, relative to Eq. (25 ), we have an additional term of 
(2£i + l)" 1 / 2 to keep the variance in each i\ fixed since 
we are now distributing power over all mi-modes. We 
adopt the fiducial redshift-dependent faint-end slope of 
the luminosity function of 



s(z) = 0.094 + 0.155z + 0.1652 2 



(32) 



estimated from the simulations of [23 , assuming a DES i- 
band magnitude limit of 24. This functional form roughly 
describes the trend that the higher redshifts are most 
affected by variations in the survey depth. We emphasize 
that this form for s(z) is meant purely for illustration, 
as different galaxy samples will have different s(z). We 
consider biases in the non-Gaussianity parameter /nl, 
the (constant) equation of state of dark energy u>, and 
the square root of the dark energy figure-of- merit, which 
is the inverse area of the 95% contour in the wo-w a plane 
[20] • Note that the latter quantity takes into account 
the temporal variation of DE, and the square root serves 
to compare it fairly to the bias in constant w; the two 
quantities, a(w) and FoM 1 / 2 , show very similar behavior 
in these results. The error bar at each l\ shows the rms 
dispersion of the (2^i + 1) values of m\ into which we put 
the systematics. So, for example, at £i=6 and for cither 
one of the rms values for the calibration error, the error 
bars show the dispersion in the bias/error ratios for 13 
different values of mi. 

Figure [3] clearly indicates that systematic errors have 
the largest impact at largest angular scales assuming a 



fixed contribution to the variance from each multipole. 
Bias in the non-Gaussianity parameter /nl is larger than 
that for the dark energy parameters, which is expected 
because most of the information on /nl comes from large 
angular scales which are particularly susceptible to cali- 
bration errors. For the calibration systematics at smaller 
scales - £\ beyond six or so corresponding to variation at 
scales less than about 30 degrees on the sky - the effect 
of the systematics asymptotes to a smaller value. Since 
calibration errors at large angular scales are the most 
damaging, it is sufficient to consider only those, and our 
choice of the maximum multipole of where the error en- 
ters, li < ^caiib,max = 20 is therefore sufficient. 

Figure [4] is similar to Fig. [3] except now we show the 
effects when the magnitude error is split equally among 
all multipoles less or equal to ^i jmax (instead of all of it 
being lumped in a single multipole as in Fig. [3]). Given 
that the effect of the systematic error decreases with l\ , 
the biases in the cosmological parameters are larger than 
in the previous figure. Qualitatively, the two figures paint 
a consistent picture of the potentially deleterious effects 
of the calibration variations even at a level corresponding 
to 0(0.001-0.01) magnitudes. 



C. Example I: Corrections to dust maps 

We now study a specific scenario of calibration system- 
atics: corrections to the SFD j3U] dust extinction maps. 
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Dust in our Galaxy causes extinction, which in turn alters 
the observed galaxy fluctuations across the sky. While 
the Galactic dust has been accurately mapped out, wc 
ask how accurately it needs to be mapped out in order 
not to bias the cosmological parameters. 

To start, we need a model for the variations in the 
SFD map. We adopt results from the work of Peek and 
Graves [37] (hereafter PG10) who used 'standard crayons' 
- objects of known color - in the SDSS to correct the SFD 
maps over the north galactic cap region (for a related 
work, see Schlafly and Finkbeiner Schematically, 
therefore 

calibration variations = (PG10 - SFD). 

The SFD map is shown in the top left panel of Fig [5] 
while the PG10 correction is displayed in the top right 
panel. To convert this Eb-v map to SN/N fluctuations 



(see Eq. ( 30 ) ) , we assume DES observations in the i-band, 
for which R = 1.595 [29]. 

The bottom left panel in Fig. [5] shows the angular 
power spectrum extracted using Polspice software pack- 
age [4] . As explained in Appendix C] at length, the most 
reliable way of modeling the calibration errors was to 
first extract the power from the map, then generate a 
full-sky realization consistent with that power (using the 
isynfast routine in HEALPix). We explicitly verified 
the intuitive expectation that the results do not depend 
much on the realization. 

The bottom right of Fig. [5] shows the resulting biases 
in /nl and the square root of the dark energy FoM, as 
a function of the faint-end slope of the luminosity func- 
tion s(z), assuming calibration variations are given by 
the PG10 corrections, and that s(z) is constant in red- 
shift. As mentioned around Eq. (28), the biases are very 
sensitive to the s(z), scaling very nearly as s(z) 2 . This 
can be easily understood: in Sec. |III C| we mentioned 
that the bias in the power spectrum is dominated by the 
added calibration power cx |Q m | 2 , while the (real-space) 
calibration field is linear in s(z) (Eq. [29]) ; hence 



(33) 



5p a oc 6Ce oc \ci m \ 2 oc s(zY 



We use a map (Jim Annis, private communication) sim- 
ulating observations over 525 night of observation spread 
over 5-years; see the top panel of Fig. [6] The observ- 
ing conditions on the site are based on historical atmo- 
spheric data of the CTIO site between 2005 and 2010. 
The survey tiling strategy uses multiple massive overlaps 
to generate as homogeneous as possible a survey. Each 
part of the sky is imaged ten times in each of the five 
DES filters (grizY). For simplicity, we only focus on the 
i-band survey-depth map. 

The effect of the unaccounted-for variability in the sur- 
vey depth is the same as that of the photometric calibra- 
tion error. However, the variability is large and expected 
to be taken into account; therefore, we (arbitrarily) adopt 
the final calibration error to be equal to one-tenth of the 
depth-variation map (i.e. 1/10 of its amplitude shown in 
the top panel of Fig. |6|. In other words, we assume 

calibration variations = — x (i-band variability map). 

We follow the same procedure as with the dust ex- 
ample above, and calculate the power spectrum of the 
depth variability map using Polspice; see the bottom 
left panel of Fig. [6] The variability of the survey depth 
will of course be accounted for in the data analysis - if 
it were not, it would lead to large biases in cosmologi- 
cal parameter estimates (as we easily verified using our 
formalism). The question is, then, to what accuracy do 
these variations need to be understood? 

We answer that question by plotting, in the bottom 
right panel of Fig. [6j the bias in the (square root of the) 
DE FoM, and non-Gaussianity /nl, as a function of the 
faint-end slope of the luminosity function s(z). As in 
the previous example of the corrections to the SFD dust 
maps, we find that the biases in the cosmological parame- 
ters are significant, and that they strongly depend on the 
faint-end slope of the luminosity function. In fact, even 
assuming that only 10% of the variability in the survey 
depth is the "calibration error" - the case shown in the 
bottom right panel of the Figure - the bias/error ratios 
are still large if s(z) > 0{1). 



Therefore, the faint-end slope of the luminosity function 
is a key factor relating the photometric magnitude vari- 
ations to the cosmological parameter biases. Steep faint- 
end slopes will lead to particularly stringent requirements 
on our understanding of the large-angle photometric vari- 
ations in the survey. 

Regardless of the value of the faint-end slope, however, 
the bottom right panel of Fig. [5] shows that the effects 
of the imperfectly estimated Galactic dust on the cosmo- 
logical parameters can be very significant. 



D. Example II: variability of survey depth 

Our second example is based on the expectations of 
photometric depth variations of the Dark Energy survey. 



VI. CONCLUSIONS 

In this paper we made a first fully general study of the 
effect of the photometric calibration variations on the 
measured galaxy clustering angular power spectra. We 
derived a general formula for how a calibration variation 
with arbitrary spatial dependence affects the measured 
galaxy angular power spectrum. We illustrated the re- 
sults assuming the standard set of cosmological param- 
eters (including /nl), DES-type dataset with five tomo- 
graphic bins out to z ma x = 1, and two specific examples 
of real-world photometric redshift calibrations. We now 
summarize our findings. 

Photometric variations modulate the observed angu- 
lar distribution of galaxy counts according to Eq. pj). 
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SFD extinction map Peek-Graves 2010 corrections to SFD extinction map 
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FIG. 5: Top left: Schlegel et al. |30| SFD extinction map, Eb-v(H), in Galactic coordinates with the 10 degrees Galactic plane 
cut. Top right: corrections to the SFD map from the work of Peek and Graves [27] ■ Bottom left: angular power spectrum 
of the PG10 map extracted by Polspice and shown without the usual 1(1 + l)/(2n) term so that the relative contribution of 
different multipoles can be more easily seen. Bottom right: bias/error ratios for /nl and the square root of the DETF FoM 
assuming PG10 map represents the calibration error, as a function of the faint-end slope of the luminosity function s. Note that 
the biases increase very sharply with s, roughly scaling as s 2 . The desired bias/error limit (horizontal dashed line) is exceeded 
already for s ~ 0.3 for /nl and s ~ 0.8 for the dark energy equation of state. 



For Peek-Graves 2010 correction to dust map 




This modulation translates into additive and multiplica- 
tive changes to the observed density fluctuation field, cf. 
Eqs. (|6| and ([9|, which in turn generate additive and 
multiplicative changes to the observed power spectrum. 

As shown in Eq. ( fTTj ), photometric variations across 
the survey masquerade as apparent violations of statisti- 
cal isotropy. Hence, explicit tests of statistical isotropy 
could provide a useful way to identify unaccounted-for 
variations in the photometry. In this paper, we focused 
on the effects in the angle-averaged power-spectrum, cf. 
Eq. (14 1. We found that large-angle modulations of 
power (dipole, quadrupole, etc), are particularly dam- 
aging to cosmological analysis. We demonstrate this ex- 
plicitly (cf. Eq. ( 26 1 and Fig. [2| for the case where the 
variance in the photometric calibration error field is con- 
centrated in one multipole £± at a time. Note that the 
spatially uniform photometric decrement or increment 
across the sky (i.e. the monopole, t\ = 0) is unobservable 
since it only affects the mean number of galaxies in the 
survey. 



Specializing in the angle-averaged power spectrum as 
done in Eq. (14), one can explicitly show that largest- 



angle fluctuations are dominant (for a fixed induced vari- 
ance on the calibration error field c(n)). (Fig. [2]). More- 
over, highest-redshift clustering measurements are most 
susceptible to the photometric variations, essentially be- 
cause their angular power is the smallest and thus is more 
affected by the photometric variation. 

Less obviously, we find that the additive errors (e.g. 



term proportional to \cg m \ in Eq. (14)) are typically 



dominant over the multiplicative biases (terms propor- 
tional to the coefficients U) for all redshift bins and at 
large angular scales. The reason is simple: because they 
couple different multipoles, multiplicative terms are sup- 
pressed relative to the additive ones by the fiducial an- 
gular power spectrum Cg factor; see the term with Cg 2 in 
Eq. (14). Since Cf < 1 even at low-z (and all £), the ad- 
ditive terms dominate the error budget if all £ modes are 
used in the analysis. However, at slightly smaller angular 
scales (£ > 10) the multiplicative error terms dominate 



13 



DES magnitude limit 




3x10 



2x10" 



u 



1x10 




10 20 30 40 50 

Multipole 1 




J 0.1 0.3 1 3 

Faint-end slope of luminosity func, s(z) 



FIG. 6: Top panel: i-band magnitude limits estimated for the upcoming observations of the Dark Energy Camera at CTIO 
as a function of angular position. The pattern of variations in the magnitude limits are set by the variations in the observing 
conditions and the survey tiling strategy over the five years of the survey. Bottom left: power spectrum of the map on the 
left, extracted using Polspice and shown without the usual 1(1 + l)/(27r) term so that the relative contribution of different 
multipoles can be more easily seen. Bottom right: biases in the cosmological parameters vs. the faint-end slope of the luminosity 
function s(z) assuming calibration error maps is consistent with a fixed fraction of 10% of amplitude (or 1% of power) of the 
magnitude-limit map shown in the top (bottom left) panel. The desired bias/error limit (horizontal dashed line) is exceeded 
for s(z) ~ 1. 



the error budget and can significantly bias the cosmolog- 
ical constraints, as discussed in Sec. |Vj Therefore, it is 
important to include both multiplicative and additive as- 
pects of the calibration error to accurately model biases 
in cosmological analyses. 

The photometric variation calibration errors affect the 
cosmological parameter measurements at largest observ- 
able scales biasing dark energy and, especially, primor- 
dial non-Gaussianity, the latter of which is intrinsically 
has signal at large scales and is particularly susceptible 
to these systematic errors. Figures [3] and |4j the principal 
plots in this paper, show these cosmological parameter 
biases for a fixed contribution of the calibration error at 
each multipole separately and for a range of multipoles, 
respectively. In Sec. |V]we further give two specific real- 
world examples of what the photometric variations could 
occur: errors in mapping the dust in our Galaxy, and 
variations in survey depth. We find that these calibra- 
tion errors lead to potentially large cosmological biases, 
especially if the faint-end slope of the luminosity function 



s = dlog iV/dm| mmax is steep. In particular, in the Fisher 
matrix approximation, the cosmological biases scale as 
s 2 . 

As a by-product of this work, we developed a reason- 
ably fast algorithm, and provide a code, to calculate the 
full biased angular power spectrum of galaxies given an 
arbitrary photometric calibration variation map on the 
sky. As mentioned above, this calculation wouldbe rather 
trivial if we could assume that the power spectrum errors 
are purely additive; however we demonstrated that multi- 
plicative errors are important and may in fact dominate, 
necessitating the more numerically intensive calculation 
for surveys of interest. 

From our analyses, it appears that the total rms of the 
calibration error has to be kept at the level at somewhere 
between 0.001 and 0.01 magnitudes rms, depending on 
exactly to how large a scale one wants to push in or- 
der to maximize extraction of cosmological information, 
in order not to bias the cosmological parameters appre- 
ciably. This is a very stringent requirement! Achieving 
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it, however, may not be as difficult as it sounds given 
that this is the time-averaged error to be tolerated at 
the end of the whole survey. Moreover, there are sev- 
eral other tools that we have at our disposal that we 
did not consider in this preliminary work. For example, 
one could use the survey itself to internally determine 
("self-calibrate") the photometric- variation errors, much 
like similar weak lensing, cluster count, or type la su- 
pernova surveys are doing or planning to do with their 
systematic errors. One could also use measurements of 
the higher-point correlation functions to help determine 
these nuisance parameters (in our language, the Q imi ). 
We leave these promising avenues for future study. 
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Appendix A: Derivation of the calibration-induced 
bias on the measured angular power spectrum 

Consider a survey where fixed, deterministic calibra- 
tion error, c(n), biases galaxy counts. In other words, 
given the true galaxy number counts in some direction 
N(h), the observed number is 



JV obs (n) = [l + c(h)]N(h). 



(Al) 



Let us expand the various quantities in harmonic space. 
First of all, true (i.e. no-systematic-errors) fluctuations 
in the galaxy counts can be expanded in the spherical 
harmonic basis 



N(n) - N 
N 



E a tm Y lm (h), (A2) 

1=1 rn=-l 



where, notice, the monopole is zero by definition of the 
left-hand side. Meanwhile, we can expand the calibration 
field relative to its fiducial value of zero (corresponding 
to no error) 



c(n) 



oo t 

EE 

e=o m=-e 



ci m Y im (h) 



(A3) 



Notice the statistical properties of the two fields 

(ai m ) = 0; (ai m a}, m t) = 8 mm iSwCe (A4) 

(q ot ) = ce m ; (c lm c* tm ,) = cg m c* e , m , (A5) 



in other words, N(h) is the Gaussian random, isotropic 
field as predicted by inflation, while c(n) is a determin- 
istic function given by calibration errors in a survey. 

Let us now calculate the observed density contrast of 
galaxies: 



<5 obs (n) 



-^obs ^obs 

JV(l + c(n)) 
N(l + e) 



- 1 



1 



1 + e 



[8(h) (1 + c(n)) + c(n) - e] 



1 + e 



8(h) ( 1 + E cimYi™ + E c t™ Y <- 



where Yg m = Y^ m (n), and we expanded the photometric 
calibration variation field c(n) into spherical harmonics. 
Here e is the relative bias in the measured mean of galax- 
ies: 



Nobs = (A r obs(n)) sk y 

= N + ^c im (N(h)Y tm (h)) sky 

= N(l + e) 
(the () s ky denotes sky average), so that 

i Y c hn(N(h)Y em (h)) 



(A7) 



sky 



(A8) 



It is now rather straightforward to calculate the mul- 
tipole coefficients of the observed density field. Writing 



5 obs (n) = ^W^ m (n) 



(A9) 



equating to the last line of Eq. ( |A6[ ) and inverting by mul- 
tiplying with Yl m and using the orthogonality relation, 
we have 



1 



a hn + Qtn + ^ ' -^mim 2 m c fimi a £ 2 ni2 

(A10) 



- V 47T£ S m S m0 



where to obtain the last term in the last line we used 
1 = Vi^Y 00 . 

Taking the two-point correlation of te m is now straight- 
forward, and things are simplified because all terms pro- 
portional to a single power of ai m (or its conjugate) van- 
ish - remember that ci m are just some numbers here. 
Moreover, we can ignore the term proportional to S(,o 
since it only affects the monopole which is not used in 



cosmological constraints. Thus we recover Eq. (11): 
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(tlmt*i> m i) — 



(1 



Jmm'OU- 



U m > m Ci> + (u mm ,y c e + 



E 



u, 



1 2 e 



(^m,m')* C«2 + c imC* e , m , 



(All) 



What about e? It can be evaluated directly in real 
space, given the calibration error map. Or it can be cal- 
culated in harmonic space, starting from its definition in 
Eq. ( A8 1 combined with Eq. (JlJ , and using the orthogo- 
nality theorem for spherical harmonics 



whose form we take to be 



(Y em Y { 



I'm' / sky 



(-i) ? 

47T 



■5u'S m {-m')- (A12) 



The result is 



coo \ ~> c-i m a lm 



'Air 



E 



■in 



(A13) 



= a lm . In cases 



where we used the identity (— l) m a^(_ r 
where / s k y < 1, the orthogonality relation in Eq. (A12) 
does not hold, but Eq. (A13) still does if the coefficients 
ai m are interpreted as the cut-sky harmonics of the den- 
sity field. 



Appendix B: Angular power spectra and their 
information content 



We model model the angular power spectra of galaxy 
density fluctuations as follows. In Limber's approxima- 
tion, the angular power spectrum is given by 



P, 



dzr(z)H(z)W t (z)W 3 (z)A 2 k 



, z 



I 

r[z) 
(Bl) 

where A 2 (fc) = k 3 P(k)/(2ir 2 ) is the dimensionless power 
spectrum, r(z) is the comoving angular diameter distance 
and H(z) is the Hubble parameter. The weights Wi are 
given by 



Wi(z) 



n{z) 



[H(z^)-H(z^ n )}, (B2) 



n(z)dz 



(i) 



were H(x) is the Heaviside step function and z^ in and 
are the lower and upper bound of the ith redshift bin 
and n(z) is the normalized radial distribution of galaxies 8 



„(*) 



8 Note that a sometimes-used alternative definition of n(z) refers 
to the spatial density of galaxies (e.g. Hu and Jain 1190 : it is 
related to the quantity we use via dN/dz = n(z) f2 r 2 (z) / H (z) , 
where Q is the solid angle spanned by the survey, and r and H 
are the comoving distance and Hubble parameter, respectively. 
Note also that our W(z) is equivalent to the quantity f(z) from 
Ho et al. [16] . 



n(z) = ^z 2 e- z ' z \ 
zz 



(B3) 



which peaks at z max = 2zq. In this work we use zq = 0.3 
so that the radial distribution peaks at redshift 0.6. 

At the lowest multipoles the Limber approximation is 
not accurate any more; see Fig. [I] At £ < 30 we adopt 
the full expression for the power spectrum; using notation 
from e.g. [H] this is 



dhikA 2 (k,z = 0)Ii(k)Ij(k) 
D(z) 



(B4) 



h{k) = b t (k) dzWi{z)^Mkx{z)) 



where x(r) is the radial distance, and x( z ) = r ( z ) m a 
flat universe that we consider. Here D(z) is the linear 
growth function of density fluctuations, so that S(z) — 
[D(z)/D(0)]S(0). At multipoles i > 30 we continue to 
use the Limber approximation since the computation is 
much faster. The power spectra are shown in Fig. [TJ Note 
that, when using Eq. (B4|, the non- linear corrections to 



clustering are negligible since we only apply this equation 
at very large scales. 

Figure [7] shows the contribution to the signal-to-noise 
squared for each cosmological parameter separately per 
multipole - in other words, we show contribution to the 
Fisher matrix element F aa for each parameter p a from 
a single multipole t. There is a clear trend of increased 
information at high multipoles, except for /nl whose in- 
formation largely comes at the lowest multipoles, despite 
cosmic variance [7]. 



Appendix C: Extraction of full-sky calibration from 
cut-sky data 

In order to estimate the effect of photometric cali- 
bration errors, we perform all calculations in multipole 
space. In particular, our formalism requires the multi- 
pole coefficients of the calibration variation field, Q imi , 
as input. In practice, this job will be left to observers 
who will have a number of tools at their disposal to esti- 
mate the photometric variation in their survey (in mag- 
nitudes, say, as a function of sky position) in their survey. 
They, or their theory friends, can then use our formalism, 
Eq. ( 14 ) , to estimate the effect of calibration variations 



on the angular power spectra and the cosmological pa- 
rameters. 
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FIG. 7: Signal-to-noise squared contribution (i.e. contribu- 
tion to diagonal Fisher matrix element with no systematics 
present) per single multipole I for each cosmological param- 
eter. Note that we use the Limber approximation only at 
t > 30 in order to more accurately capture information on 
/nl at £ < 30. The reason for the small discontinuity at 
I = 30 is due to the switch from full to Limber power spec- 
trum. 



However it is nontrivial to obtain the coefficients ci imi 
given the cut-sky observations, i.e. partial sky coverage: 
recall, for example, that the Peek-Graves dust-correction 



map from Sec. VC covers about 1/4 of the sky, while 
the DES depth map from Sec. |V D| covers only 1/8 of 
the sky. Note that the naive reconstruction from the 



observed map, Q rmi = J Y e * imi (ii)SN/N(h)dQ, would 

give the cut-sky multipoles which describe the field that 
is zero outside of the observed region, which we do not 
want: the cut-sky multipoles have wrong amplitudes and 
harmonic structure relative to the full-sky 'truth'. 

Filling in the missing sky is a well-known problem in 
the cosmic microwave background literature. While the 
reconstruction of the power Cg is relatively straightfor- 
ward, reconstruction of the temperature field ag m is chal- 
lenging and typically works well only when substantial 
portions of the sky (/ 8 ky ^ 0.8) are observed. Fortu- 
nately, in this work we are only interested in obtaining 
approximate values for the ci imi , consistent with the cut- 
sky map of the calibration field, in order to estimate its 
effects on cosmological parameters. 

We have attempted, without success, two rather well- 
known techniques for reconstructing the true, full-sky 
cg imi from partial sky observations: 

• We tried the direct reconstruction of the field, us- 
ing the maximum-likelihood procedure (see e.g. 
O [lOl [E]); this approach is used, for example, 
to reconstruct the CMB power at i < 30 [25] . 
While this approach returns the true structure in 
the observed area of the sky (as it must), the recon- 



structions add spurious structure and extra power 
- clearly visible by eye - in the unobserved areas 
of the sky. This remains the case even after exten- 
sive testing and fine-tuning of the procedure, and it 
only returned trustworthy results at high sky cov- 
erage (/ s ky > 0.8). Therefore, this approach fails to 
extrapolate the photometric-variation field beyond 
the observed region even approximately. 



• We tried "harmonic inpainting" [181 EI] of our two 
maps, were the missing portion of the sky is filled 
in assuming that the density field is Gaussian, and 
starting with a guess for the covariance matrix of 
the field (which is just given in terms of the theo- 
retical angular power spectrum). As with the di- 
rect reconstruction, we find large spurious power in 
the unobserved portion of the sky in realistic cases 
when sky coverage is small, / s k y <C 1. The negative 
results persist even after attempts to fine-tune the 
analysis to the maps at hand, again in parallel with 
the reconstruction method. 

Given these failures and the fact that we only need a 
rough estimate of the full-sky calibration field cg imi , we 
resort to the simpler scheme of reconstructing the large- 
angle power spectrum and drawing realizations of the 
field consistent with it. In particular, we 

1. Calculate the power spectrum using Polspice 
package [J] to calculate the angular power spectrum 
from the cut-sky map. Polspice uses the pixel- 
based approach to calculate the real-space angu- 
lar power first and then, with appropriate apodiza- 
tions, convert it to the Cg. 

2. Sample from these Ci t , i.e. cg l0 = Af(0, G^) and 
4^ = M{0,C t j2) for To! ^ 0, where 7V(m,<7 2 ) 
is the Gaussian normal function with mean \x and 
variance cr 2 . 

3. Repeat the previous step a number of times 
to ensure appropriate averaging over realizations, 
though in practice we find that the each realization 
(for a fixed estimated Cg t ) leads to similar results. 

The approach returned a smooth power spectrum at large 
scales (see the Cg panels in Figs. [5] and [6]), and therefore 
reliable multipole coefficients describing calibration vari- 
ations across the whole sky. 

The reader might be worried that, in the above- 
mentioned successful approach, we assumed Gaussianity 
for the calibration field by assuming that the Cg t de- 
scribes the full-sky field, and moreover assuming that 



the 



C£ l7 



are Gaussian random variates. This is defi- 



nitely an assumption needed to get the coefficients in our 
particular examples described in this Appendix but not 
in general. Moreover, we are not concerned about the 
assumption of Gaussianity of the calibration-variation 
field: Eq. (Ill, for example, illustrates that the calibra- 



tion field will break statistical isotropy of the observed 
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density fluctuations regardless of the nature of the sta- 
tistical distribution of the ensemble from which th ci imi 
come from. Modest dispersion among the different ra\- 
modes in Fig. [3] confirms that the principal property of 
these coefficients that drives their effect on cosmology are 
their amplitudes as a function of t\, not their phases. 

We emphasize here that the somewhat-challenging task 
of needing to determine the (full-sky) ci lTni discussed in 
this Appendix is largely orthogonal to the overall goal of 
this paper of quantifying the effects of photometric varia- 
tions. In practice, the observers will have an opportunity 
to estimate the photometric calibration variations across 
the sky in real space using a number of methods (e.g. 
using theoretical estimates, or measurements across part 
of the sky as we did in our examples). Converting from 
c(n) to the full-sky ci imi is then just a purely mathe- 
matical exercise in which only an approximate answer is 
required. 



while the latter quantity we now evaluate. 

At a fixed multipole of the true density field £2 and a 
calibration error at some multipole £%, only the observed 
multipoles in the range max{|4 — 4|,0} < £ < £2 + £1 
are affected; this is well known feature of the coupling of 
angular momenta in e.g. quantum mechanics. Therefore, 
for a fixed values of I and m, the total number of (£2, m.2) 
pairs that can possibly lead to nonzero U coefficients is 



naively 



E 



(2£ + 1). However, we must 



also remember of a selection rule 9 that mi + iri2 = m, 
which means that for a fixed observed-field coefficient m 
and calibration-field coefficient mi, only one coefficient 
m,2 (instead of 24 + 1 of them) survives. Therefore, for 
a fixed £, the total number of U coefficients to tabulate 



(e.g. for use in Eq. (14)), is of order 



Appendix D: Tabulation of coefficients 



In order to speed up calculations of the observed galaxy 
overdensity field tg m and its correlations in harmonic 
space, we need to calculate a very larg e n u mbe r of ge- 
ometrical quantities ]?^ 2 j 2m ; see Eqs. ( 10 )-( 12 ). Recall 
that each (£, m) pair has 2£+l quantities lor a fixed m, or 
approximately ^ ax quantities for all £ up to some £ max . 
Calculating all of the R coefficients would then naively in- 
volve 0(^ lax ) operations which, for 4iax ~ 1000, would 
be 10 18 calculations of the Wigner-3j symbols! Calculat- 
ing even a fraction of such a large number of coefficients 
is clearly unfeasible. 

We bring the computation required to a manageable 
size and speed as follows. We decide to tabulate the co- 
efficients U for the minimum required number of quadru- 
plets of its indices, and carry out the summation in 
Eq. ( 12 ) for each U. Moreover, since our goal is to con- 
centrate on large-angle (low-4) systematics, which lim- 
its the maximum multipole of the systematics, £\ (see 

Eq. (Il2|)), we Choose £1 < 4alib,max with 4alib,max = 20. 



~ (2£+l)(2£, 



calib.max 



+ l)JVi 1 ,mi 



Therefore, the total number of evaluations of the Wigner- 
3j symbols will be approximately equal to the number of 
(£1, mi) pairs times the number of the U coefficients; the 
former quantity is 



£1=0 



(24 + 1) = (Calib.max + I) 2 , (Dl) 



calib.max 



(D2) 



where the factors of 2£ + 1 and 24aiib,max + 1 refer to 
the number of m coefficients and (£2,^12) pairs, respec- 
tively. So for a typical calculation with 4aiib,max — 20, 
the number of evaluations for a single multipole £ ~ 250 
is of order 10 7 ; for of order 30 values of £, including the 
densely spaced coverage of lowest observable multipoles 
where the effects of photometric variation are the largest, 
this is still less than O(10 9 ) evaluations of the Wigner-3j 
symbols, which is entirely feasible and takes of order a 
minute on a single 16-core desktop computer. 

The last detail is actual computation of the Wigner- 
3j symbols: we use the Gnu Scientific Library (GSL) 
routines at I > 50. At higher multipoles these break 
down, and we use a special-purpose routine 10 that 
uses approximate expressions for the Wigner coefficients. 
These expressions hold provided £,£2 3> £1, but given 
that £ > 50 in the regime we are using this approxi- 
mation and we only consider l\ < 4aiib,max — 20, this 
condition always holds. 
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